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We consider the proton and neutron quasiparticle orbits around the closed-shell '^Ni and ''^Ca isotopes. It is 
found that large model spaces (beyond the capability of shell-model applications) are necessary for predicting 
the quenching of spectroscopic factors. The particle-vibration coupling is identified as the principal mechanism. 
Additional correlations — due to configuration with several particle-hole excitations — are estimated using shell- 
model calculations and generate an extra reduction which is < 4% for most quasiparticle states. The theoretical 
calculations nicely agree with {e,e'p) and inverse kinematics knockout experiments. These results open a new 
path for a microscopic understanding of the shell-model. 
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Introduction. — Single-particle states at the Fermi surface 
of shell closures (or quasiparticles) play a crucial role in 
nuclear structure. Their strengths [the spectroscopic factors 
(SFs)] and energies are essential inputs in large-scale shell- 
model (SM) calculations, which provide the most success- 
ful description of medium mass Qsld and Oplf nuclei Ol. 
Single-particle wave functions and SFs also enter the nuclear 
matrix elements that describe nucleon capture and emission in 
stellar burning processes j^]- 

The need for understanding quasiparticles has become 
more compelling in recent years due to the realization that 
their energies evolve, with changing_proton or neutron num- 
ber, due to the nuclear tensor force |3, 4]. The quenchings of 
absolute spectroscopic factors (SFs) have also been observed 
to change dramatically at the dripUnes 01 • These effects 
cause the breakdown of conventional magic numbers far from 
the valley of stability and can give rise to new exotic modes 
of excitation. In order to predict the properties of nuclei at the 
limits of stability, it is imperative to derive robust models for 
quasiparticle properties. Recently it has been suggested that 
relative SFs among small quasiparticle fragments can be de- 
scribed within the SM approach ff\. This Letter discusses the 
problem of the absolute quenching of SFs and shows that this 
requires accounting for coupling to collective resonances and 
large model spaces, beyond the capability of the SM. 

The reactions for transfer of a nucleon to/from the initial 
state I'Pg ) depend on the overlap wave function 
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where a can label either particle or hole states. SFs are iden- 
tified with the normalization integral of iA^'^Hr) and give a 
"measure" of what fraction of the final wave function, l*!*^^'), 
can be factorized into a (correlated) core plus an independent 
particle or hole. Strong deviations from the independent par- 
ticle model (IPM) — that is, a Slater determinant with fully oc- 
cupied orbits — signal substantial correlations and imply the 
onset of non trivial many-body dynamics. For stable nuclei, 
a large body of data has been accumulated from {e,e'p) ex- 
periments. These studies showed that proton absolute SFs for 
isotopes all across the nuclear chart are uniformly quenched 
to 60-70% of the IPM value E 1 in . Short-range correlations 



(SRC) — originating from the repulsion of the nuclear force at 
short distances — have been investigated as a possible cause. 
However, they explain only a small fraction of the observed 
reduction fill EH], as confirmed by recent experiments at Jef- 
ferson Laboratory iflA 15]. Thus, one is left to understand the 
role long-range correlations (LRC) eff'ects (see, for example, 
Refs. lilEl and references therein). 

Successful ab-initio calculations have been reported for 
light nuclei using the variational Monte Carlo method lIlSll . 
For larger masses, self-consistent (SC) Green's function 
(GF) theory |.16. 19] gave encouraging results for '^'^Ca and 
^"Zr ll20[l . However, these early calculations focused only on 
one particular type of correlations [either coupling to particle- 
hole (ph) or particle-particle and hole-hole phonons]. Re- 
cently, we have improved the SCGF method 12111 to include 
correlations from the remaining collective modes, extended 
it to large-scale no-core calculations up ^^Ni, and based our 
studies on realistic nuclear forces |^, 23). Thus, we are now 
in the position of making accurate studies of the eff'ects of cor- 
relations. In addition, the SCGF equations describe directly 
the interactions among quasiparticle states, which allows ex- 
tracting single-particle orbits and effective interactions for the 
SM. This Letter applies these new capabilities to study in de- 
tail quasiparticle states in ^^Ni and "^^Ca and to extract an ef- 
fective Hamiltonian for Qsld and Oplf orbits in the SM. The 
SCGF results will then be compared to the SM ones to in- 
vestigate the role of all dififerent types of correlations on the 
quenching of SF. 

Formalism. — The calculations were performed in a large 
model space including of up to 10 major harmonic oscillator 
shells. All details are reported in Ref. [23:] and here we give 
an overview of the points needed for the present discussion. 
Spectroscopic factors were extracted from the one-body prop- 
agator (spin/isospin indices are omitted for simplicity) il9ll . 
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where the residues are the overlap amplitudes ([TJ and the 
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FIG. 1: (Color online) Left. One of the diagrams included in the correlated self-energy, Arrows up (down) refer to quasiparticle 

(quasihole) states, the n*''''' propagators include collective ph and charge-exchange resonances, and the g" include pairing between two 
particles or two holes. The FRPA method sums analogous diagrams, with any numbers of phonons, to all orders i2ll. [2^ . Right. Single- 
particle spectral distribution for neutrons in ^*Ni, obtained from FRPA. Energies above (below) Ef are for transitions to excited states of 
^'Ni ('^Ni). The quasiparticle states close to the Fermi surface are clearly visible. Integrating over r [Eq. l[4j] gives the SFs reported in Tab.|l] 



poles give the experimental energy transfer for nucleon pickup 
(knockout) to the excited states of the systems with Ah-1 (A-1) 
particles. The propagator (|2]i is obtained by solving the Dyson 
equation {gicS} - ^'***(w) -i- g^^\(jL)) 2*((y) gico)], where 
g*"'((y) propagates a free nucleon. The information on nuclear 
structure is included in the irreducible self-energy, which was 
split into two contributions: 



E*(r, r'; u) = S*"'(r, r'; u) + £(r, r'; u) 



(3) 



The term YJ^^io)) includes both the nuclear mean field (MF) 
and diagrams describing two-particle scattering outside the 
model space, generated using a G-matrix resummation 1241. 
As a consequence, it acquires an energy dependence which 
is induced by SRC among nucleons |23]. The second term, 
£(tij), includes the LRC. In the present work, £(w) is calcu- 
lated in the so-called Faddeev random phase approximation 
(FRRA) of Refs. ll2lll25ll . This includes diagrams for particle- 
vibration coupling at all orders and with all possible vibration 
modes, see Fig.[Tl as well as low-energy 2plh/2hlp configu- 
rations. Particle-vibration couplings play an important role in 
compressing the single-particle spectrum at the Fermi energy 
to its experimental density. However, a complete configura- 
tion mixing of states around the Fermi surface is still missing 
and would require SM calculations. 

Each spectroscopic amplitude '(r) appearing in Eq. ^ 
has to be normalized to its respective SF as 
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where S*.(w) =< 0'o,|S*(w)|i^„ > is the matrix element of 
the self-energy calculated for the overlap function itself but 
normalized to unity (J dr |i^Q,(r)p = 1). By inserting Eq. (O 
into (01), one distinguishes two contributions to the quenching 
of SFs. For model spaces sufficiently large, all low-energy 
physics is described by £(<y). Then, the derivative of ll'^'^iiij) 



accounts for the coupling to states outside the model space 
and estimates the effects of SRC alone 13311 . 

In general, the SC self-energy (|3]l is a functional of the one- 
body propagator itself, 2* = Hence the FRPA equa- 
tions for the self-energy and the Dyson equation have to be 
solved iteratively. The mean-field part, ll^^lg}, was calcu- 
lated exactly in terms of the fully fragmented propagator (|2]l. 
For the FRPA, this procedure was simplified by employing the 
Y.{g"''^~\ obtained in terms of a MF-Iike propagator 
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which is updated at each iteration to approximate Eq. (|2]i with 
a limited number of poles. Eq. (|5]l defines a set of undressed 
single-particle states that can be taken as a basis for SM ap- 
plications. This feature will be used below to estimate the im- 
portance of configuration mixing effects on the quenching of 
spectroscopic factors. The present calculations employed the 
N3LO interaction from chiral perturbation theory |26] with a 
modification of the tensor monopoles to correct for missing 
three-nucleon interactions I27ll . 

Results. — The calculated single-particle spectral function 
{S^cj^iir, Lo) - ^l^(r - r'; is shown in Fig.[T]for the case 
of neutron transfer on ^^Ni. This picture puts in evidence the 
quasiparticle and quasihole states associated with valence or- 
bits in the Qplf shell. The corresponding SFs are reported 
in Tab. U including both protons and neutrons. The first col- 
umn is obtained by including only the derivative of Yl*^^{io) 
when calculating Eq. (|4|i. Since N3LO is rather soft com- 
pared to other realistic interactions the effect of SRC is rela- 
tively small. From other models one could expect a quenching 
up to about 10% 1 16], as confirmed by recent electron scatter- 
ing experiments \ \A, ilSt |28[] . This difference would not affect 
sensibly the conclusions below. The complete FRPA result for 
SFs is given in the second column. For the transition between 
the ^''Ni and ^^Ni ground states, our result agrees with knock- 
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out reactions with fast beams within the given error bar ll29ll . 
Note that ^^Ni is close to the stabiUty Une, where the analysis 
of these reactions has been tested against electron scattering 
results. The remaining orbits have similar quenchings and are 
in line with results expected for stable nuclei 1 10]. This shows 
that the dominant mechanism in reducing the spectroscopic 
strength is to be looked for in LRC involving configurations 
among several shells around the Fermi surface. 

We now turn to the question of what are the relevant low- 
energy degrees of freedom and the importance of config- 
urations with multiple ph excitations. Excitations beyond 
2plh and 2hlp are partially included by the SC-FRPA, from 
which one can expect an accurate description of quasiparti- 
cles around shell closures. This was checked by comparing 
SC-FRPA and large-scale SM calculations in the 1 pQf shell. 
The basis and the interaction were constructed using the "un- 
dressed" single-particle orbits ^^(r), defined by Eq. (|5]l, and 
the same G-matrix interaction used in the full-space FRPA 
calculations. Shell-model calculations used the ANTOINE 
code fll,'30^ and included up to 6p6h excitations for the A=57 
systems and up to 8p8h for A=55. The resulting SFs are 
shown in the last three columns of Tab. U As an example, the 
quenching for the 793/2 ground state of ''^Ni, is predicted to be 
0.82 by FRPA and 0.79 by the SM. Comparing with the full- 
space result of 0.65, one infers that about half of the quench- 
ing is driven by degrees of freedom from outside the 1 pQf 
space. In particular, giant resonances originating from Ifia) 
and Ihco excitations carry enough collectivity to influence the 
quenching of orbits close to the Fermi surface. The difference 
between SM and FRPA in the Opl f space is AZ„=0.79-0.82=- 
0.03, which is smaller than the typical model dependence of 
SFs extracted from experiments. One can treat the AZq,s as 
small perturbations and add them to the full-space FRPA re- 
sults, leading to our final predictions for absolute SFs (third 
column of Tab.|l]l- Relative SFs with respect the other states 
in ^^Ni have been studied by re-analyzing the world {d,p) 
data fT]. These results agree with the SM and, hence, with 
FRPA ones (see Tab.|I]l. However, the model dependence of 
this analysis appears too large to constrain absolute values of 
SFs. 

Table reports the results of the same analysis for ^'^Ca. 
Also in this case, degrees of freedom from outside the first 
10 oscillator shells have a small effect, of 3-6%, on the abso- 
lute SFs. The SC-FRPA coupling to collective modes reduces 
these to about 80% of the IPM, and even lower for proton 
holes to ^^K. The SM calculations employed the "undressed" 
single-particle orbits (|5]l generated by the large space SC- 
FRPA for this nucleus. We adopted a model space (which will 
be named sd^pf) consisting of the (Lsi/2,0c/3/2, 1/53/2,0/7/2) 
orbits for protons and of (lpi/2, 1/^3/2,0/5/2,0/7/2) for neutrons. 
This can be handled by present day computers with up to 
lOplOh configurations and includes the most important exci- 
tations across the Fermi surface. The SM results in the sd^^p f 
model space do not give additional quenchings for the valence 
orbits in the 1 pOf shell. Hence, multiple ph configurations 
are suppressed even more than for the ^^Ni case. Note that 
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TABLE I; Spectroscopic factors (given as a fraction of the IPM) for 
valence orbits around ^'Ni. For the SC-FRPA calculation in the large 
harmonic oscillator space, the values shown are obtained by includ- 
ing only SRC, SRC and LRC from particle-vibration couplings (full 
FRPA), and by SRC, particle-vibration couplings and extra corre- 
lations due to configuration mixing (FRPA+AZ„). The last three 
columns give the results of SC-FRPA and SM in the restricted 1 pQf 
model space. The AZ^ are the differences between the last two re- 
sults and are taken as corrections for the SM correlations that are not 
already included in the FRPA formalism. 

the vO/5/2 quasiparticle in "^^Ca splits in two main fragments 
but their total strength remains unchanged. Experimental SFs 
for the ground (yQp-iji) and first excited (vO/?i/2) states, have 
been extracted from a combined analysis of {d,p) and {n,y) 
reactions jjlll . 

Proton-hole states in "*^K are a case in which SM correla- 
tions start becoming important. The FRPA in the sdip f space 
describes about a half of the quenching of the 0(^3/2 orbits 
in the full space. The SM calculations induce a further re- 
duction of si8%. Adding this estimate brings the 0ii3/2 SF 
to 0.59, in agreement within the experimental {e,e'p) value 
of 0.56(5) il III . Shell-model correlations are even more sub- 
stantial for the lsi/2 fragments: here the sd^pf space is suffi- 
cient to reproduce almost all of the quenching obtained in the 
largest space. Configurations not included in the FRPA are 
also more important and induce a further reduction of up to 
16% for the ground state of "^^K. The total predicted strength 
of the two lii/2 peaks is 0.50, somewhat lower that the {e,e'p) 
value of 0.61(5). 

Utsuno and collaborators have also calculated the spectral 
strength for proton removal from "*^Ca and found that the spin- 
orbit spitting between the Qd states and the ordering of single- 
particle energies are sensitive to the tensor force i32tl . These 
calculations also show that the fragmentation pattern in "*^K 
is well reproduced by the SM, as suggested in 101 ■ The ex- 
perimental strength is however lower by an overall factor of 
aiO.73. This is nicely explained by the combined SC-FRPA 
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0.56(4) 
0.53(4) 
0.08(1) 



0.87 
0.64 
0.34 



0.79 
0.48 
0.23 



-0.08 
-0.16 
-0.10 



TABLE II: Same as Tab.|l]but for ^'*Ca. The vO/5/2 particle (nlsi/2 
hole) fragments in '"Ca C"K) with lowest excitation energy are in- 
dicated by an asterisk. The experimental SFs for ""Ca are from 
Ref. those for "'K are form Ref. LUJ. 



and SM results, where the final corrected strengths of TablU 
are smaller than the pure SM ones by 0.75 for the Od3/2 peak 
and by 0.70 for the summed I51/2 strength. 

Summary. — In conclusion, we have studied the effects of 
correlations on the quenching of absolute spectroscopic fac- 
tors. The present work does not support the assumption that 
the uniform quenching observed in stable nuclei originates 
from a unique mechanism such as short-range correlations. 
SRC have been predicted [12, 16] and experimentally con- 
firmed 1 14 , Tstl to influence just a small fraction of the to- 
tal strength. We find that spectroscopic factors for valence 
orbits result from a combination of configuration mixing at 
the Fermi surface and coupling to collective resonances. The 
former can usually describe the relative spectroscopic factors 
(hence, the fragmentation pattern at low energy) and it will 
be crucial in open shell nuclei. On the other hand, the over- 
all quenching of absolute spectroscopic factors requires both 
the coupling to collective modes at higher energies and large 
model spaces that cannot be approached by shell-model calcu- 
lations. This work reported on the first consistent calculation 
of all these effects in the medium-mass region and used the 
same realistic interaction. 

It was also found that quasiparticles around shell closures 
can be described accurately by self-consistent Green's func- 
tion theory in the FRPA approximation. Moreover, the theory 



links the interactions among quasiparticle states with the un- 
derlying nuclear forces. This opens a new possibility for de- 
riving single-particle orbits, effective charges and interactions 
for the shell-model form realistic nuclear interactions. 
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